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Abstract 

We explore effects of the Shakura-Sunyaev a-viscosity on the dynamics and oscillations of slender 
tori. We start with a slow secular evolution of the torus. We show that the angular-momentum profile 
approach the Keplerian one on the timescale longer than a dynamical one by a factor of the order of 
1/a. Then we focus our attention on the oscillations of the torus. We discuss effects of various angular 
momentum distributions. Using a perturbation theory, we have found a rather general result that the 
high-order acoustic modes are damped by the viscosity, while the high-order inertial modes are enhanced. 
We calculate a viscous growth rates for the lowest-order modes and show that already lowest-order inertial 
mode is unstable for less steep angular momentum profiles or very close to the central gravitating object. 
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1. Introduction 

Stability and propagation of waves in viscous accretion 
disks have been studied extensively for more then thirty 
years, since the pioneering work of Kato (1978). Adopting 
the local approximation, he has shown that the horizon- 
tal p-mode oscillations of the Keplerian isothermal disk 
become over-stable when a small viscosity of the fluid is 
introduced. The instability arises because the azimuthal 
component of the viscous force varies in phase with the 
azimuthal velocity during oscillations and thus a positive 
work on the oscillations is done. The source of the energy 
for the oscillations is the angular momentum flowing ra- 
dially in the background flow. Under suitable conditions, 
given mainly by the temperature and density dependence 
of viscosity and also by the mode eigenfrequencies and the 
shape of the modal eigenfunctions, this process may pro- 
vide energy sufficient to amplify the oscillations against 
other normal viscous damping processes. The instability 
of p-modes is even stronger in the relativistic disk, where 
the gradient of the angular velocity is steeper. 

Later on, this problem has been reconsidered by other 
authors and the analysis has been expanded also to other 
types of oscillation modes. For example, Blumenthal et al. 
(1984) studied stability of p-modes in disks with different 
ratios of the gas to the radiation pressure, Kato (1991) dis- 
cussed possibility of viscous excitation of the one-armed 
corrugation mode in relativistic disks. The local approx- 
imation was released in the work of Nowak & Wagoner 
(1992), who estimated growth rates of p-modes trapped 
near the inner edge of the disk and g-modes trapped near 
the radius of the maximum of the radial epicyclic fre- 



quency. 

Some of the analysis go even beyond the standard 
Shakura & Sunyaev (1973) prescription that is based on 
the mixing-length theory. Nowak & Wagoner (1992) con- 
sidered the effects of an anisotropic viscosity and showed 
that the growth of the trapped g-modes is enhanced while 
that of the p-modes is suppressed. Perhaps even more ad- 
vanced is the analysis of c-modc oscillations in relativistic 
disks by Kato (1994) who uses the transport equations in 
calculating the stress tensor perturbations based on his 
second-order closure scheme. 

In this note we study effects of viscosity on oscillation 
modes of a slender torus. The main difference of this con- 
figuration from the thin disks studied in the work cited 
above are the angular momentum profiles that can be sig- 
nificantly different from the Keplerian distribution and 
the following substantial pressure and density gradient in 
the radial direction (i.e. in the same direction as the an- 
gular velocity is changing). The plan of the paper is as 
follows. In section 2 a basic theory of inviscid tori and 
the limit of slender tori are reminded. The section 3 deals 
with the secular evolution of the torus that arises when 
a small viscosity is introduced to the flow. The section 4 
is devoted to the general discussion of oscillation modes 
and the effects of viscosity on the oscillation frequencies. 
There we derive a general formula for the eigenfrequency 
shifts which we subsequently apply in sections 5 and 6. 
The section 5 contains results for the modes of high or- 
der and for Keplerian tori and the section 6 focuses on 
the lowest-order modes. Finally, section 7 is devoted to 
discussion and conclusions. 
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2. Inviscid slender tori 

The construction of an inviscid axisymmctric poly- 
tropic slender torus is described by Blaes (1985) in the 
Newtonian gravitational field and by Blaes et al. (2007) 
in the case of a general axisymmetric gravitational poten- 
tial. The fluid of the torus is in pure rotation, e.g. its 
velocity in the polar coordinates {r, 4>,z} can be expressed 
as v 1 = flS^, where il is the angular velocity. The shape 
of the equi-pressure and equi-density surfaces coincide be- 
cause of the polytropic equation of state, p cx p 1+1 / n , and 
can be uniformly described by the Lanc-Embden func- 
tion f(x). The density and pressure are then given by 
p = pof n (x) and p = pof n+1 {x) (the subscript zero indi- 
cates the evaluation at the center of the torus, where the 
density and pressure are maximal). Another consequence 
of using the polytropic equation of state is that the an- 
gular velocity is a function of the radial coordinate only, 
O = O(r). In addition, throughout the paper, we assume 
that this function is slowly varying, i.e. <9 r f2 ~ fi/r. The 
structure of these surfaces follows from the poloidal part 
of the Euler equation 



P " 



(1) 



where 1 = 2(n+ l)po/(/0o'"o^o) i s the slenderness param- 
eter j3 that gives the thickness of the torus, both the radial 
and vertical extend of the torus are of the order of /3ro . 

The limit /3 — > describes slender tori. The surfaces of 
constant pressure and density have elliptic cross-sections, 
whose centers correspond to the maximal-pressure circle 
at r = ro in the equatorial plane. It is convenient to intro- 
duce dimensionless coordinates x and y contracting with 
the torus as /3 — > 0, 

*=^, V=j- (2) 
pr pr 

The shapes of the equipressure surfaces are then given by 

f(x,y) = l-(u 1 2 -R 2 )x 2 -Q 2 z y 2 , (3) 

where ui r and lj z are the radial and vertical epicyclic fre- 
quencies normalized by flo and R 2 — (dln€ 2 /dlnr)o de- 
scribes distribution of the angular momentum £ = r 2 f2(r), 
which is in the slender-torus limit linear. It follows from 
equation (1) that fio = Ok(?*o)- The outermost ellipse 
given by / = describes the surface of the torus. On the 
other hand, the torus center at (x,y) = (0,0) corresponds 
to / = 1. 

3. Secular evolution of viscous tori 

Introduction of a viscosity causes a slow secular evolu- 
tion of the torus governed by the Navier-Stokes equation 

^- + v k V k rf + -V l p + V*$ = -V k a hl = J-; isc (4) 
at p p 



The stress tensor of the viscous flow is given by 

a lk = V (VV + V l v k ) + U- |^ (V • v)g ik , (5) 



where rj and £ are coefficients of dynamic and bulk vis- 
cosity, respectively, and g lk is a metric tensor. We extend 
the standard 'ap'-paramcterization of the shear viscosity 
77 traditional in the accretion theory to include also the 
bulk viscosity £ by introducing a coefficient a, 



rr- 



ap 



ap 



a< 1. 



(6) 



The velocity of the fluid is still dominated by the az- 
imuthal component, but also small poloidal components 
appear due to the viscosity. In our analysis, we take a as 
a small parameter with respect to which both the equa- 
tions and their solutions are expanded. For the poloidal 
velocity we assume that v l /(r£l) = O(a), i = r,z, the az- 
imuthal velocity is of the order of the Keplerian orbital 
velocity at the center of the torus. From a balance among 
leading-order terms in the poloidal part (i =r,z) of the 
Navier-Stokes equation we obtain again 



-y l p + V l ^-rVi 2 5 l r = 0, 
P 



(7) 



while the toroidal part and the continuity equation gives 
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(8) 



(9) 



The right-hand sides of the latter two equations are of the 
first order in a, these equations therefore describe a slow 
evolution of the angular momentum and the density on 
the time scale ~ l/(afi). On the other hand, the equa- 
tion (7) states the equilibrium among the pressure, grav- 
itational and centrifugal forces that is reached on much 
shorter timescale ~ I/O. This equation determines the 
structure of the equipressure surfaces in the torus and is 
identical to the case of the inviscid flow. The overall evo- 
lution can be therefore regarded as a sequence of inviscid 
tori whose shapes are given by Lanc-Embden function / 
defined in equation (3), but whose main determining pa- 
rameters R, j3 and ro are slowly changing in time. 

Calculating the total mass A4 and angular momentum 
C of infinitely slender torus by integrating p and p£ over 
the torus volume, we obtain the relation 



£ = £ K (r )M. 



(10) 



As both quantities has to be conserved during the process, 
we may conclude that the position of the torus center at 
r = ro stays unaffected. This is an artefact of the reflec- 
tion symmetry with respect to x — that appears in the 
limit (3 — > 0. The centers of thicker tori that do not have 
this symmetry will be, in general, slowly drifting in time 
toward the central objects. 

Substituting the approximate expressions (3) for / with 
R = R(t) into equations (8) and (9) and assuming a linear 
profile of £(r), we obtain 
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Fig. 1. Top: A slow secular evolution of a slender torus from 
initially constant angular momentum (re = 0) towards the 
Keplerian distribution (R r =Gi r ). In the calculation the value 
of the radial epicyclic frequency has been taken to correspond 
to the radius ro = lOAf in Schwarzschild spacctime. Bottom: 
Shape of the torus at various times (t X aQg = 0, 0.1, 0.2, 0.4, 
0.8 and 1.5). 
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Now, making an ansatz 
v r (t,x,y)=j3xV r {t), 



v z (t,x,y) = /3yV z (t), 



(13) 



we obtain polynomial equations in x, y. Comparing ap- 
propriate coefficients, we finally arrive at a set of four 
equations 
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+ V'](o 



[V r + {2n+l)V z ](J z 



-R 2 V r 
ro 



for unknowns R(t), (3(t), V r (t) and V z (t). The time of the 
secular evolution from the constant-angular momentum 
torus to the configuration characterized by R is 

rR 2 (t) 



ailo 



{A-k){Q 2 -kf 



an (4-u> 2 ) 2 
(l-q)R 2 i 



(Aq-d 2 )\n 
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A(uj 2 - R 2 ) 



0J~ 



(18) 



where q = (2n + 3)/(4n + 4). 

The solution R(t) is shown in the top panel of Figure 1. 
The angular momentum from the inner part of the torus 
is transported outward; its distribution becomes steeper 
approaching the Keplerian one (corresponding to R = u> r ) 
at long times. A characteristic timescale of the viscous 
diffusion is 



-if)- 1 



(19) 



The bottom panel shows the shape of the torus at sev- 
eral moments during the secular evolution. The torus be- 
comes wider in a radial direction and shrinks in the ver- 
tical, as the angular momentum distribution approaches 
the Keplerian one. 

In the next sections we calculate damping rates of the 
lowest-order oscillation modes due to viscosity in tori with 
arbitrary angular-momentum distributions. Since the os- 
cillations occur on timescales that are much shorter than 
iviso we will ignore the viscous diffusion in the rest of the 
paper. 

4. Perturbations 

The unperturbed state is axisymmetric and can be 
regarded as stationary on timescales shorter than t v isc- 
Therefore the t and ^-dependence of the perturbation is 
ocexp[— i(u)t — Tncj))], where u> is eigenfrequency of the oscil- 
lation mode and m its azimuthal wavenumber. Papaloizou 
& Pringle (1984) showed that it is convenient to express all 
perturbations in terms of a single quantity W = —Sp/ (per)., 
where a — u> — mCl. Linear Euler perturbation of the con- 
tinuity equation gives 



--^-(rf n Sv r ) + ^-(f n Sv z ) 
r or oz 



imf n Sv' t> 
2na 2 



P 2 



f"- l W = 0, (20) 



where a = u/flo- The perturbation of the velocity 8v l 
can be calculated from W using perturbed Navier-Stokes 
equations, which takes the form 
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The perturbation of the viscous force in terms of the ve- 
locity perturbation W is 
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where we keep only leading-order terms in the limit /3 — > 
(i.e. we neglect effects of the azimuthal curvature of the 
torus). The equations (21)-(23) together with (24)-(26) 
are then solved for poloidal velocity in terms of W. This 
procedure is easier if we assume that dW/dx > 1, i.e. that 
the wavelength of a perturbation is at least comparable 
to the size of the torus. Then, up to the first order in a, 
and in the lowest order in /3, the solutions are 

6v r 
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and 
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Substituting equations (27) and (28) into the continuity 
equation (20) and keeping only dominant terms when (3 
0, we obtain a single operator equation for W, 
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(32) 



When a = 0, equation (29) is the second-order eigen- 
value problem and describes oscillation modes of the in- 
viscid slender torus. If, in addition, R = (constant angu- 
lar momentum tori), the equation (29) is reduced to the 
first-order eigenvalue problem for the operator B. In that 
case the operator B is simply related to the operator L^°> 
introduced in Blaes et al. (2006) by B = 2nU°\ Blaes 
(1985) and Blaes et al. (2006, 2007) showed that this op- 
erator is Hermitian with respect to the scalar product 



(U,V)= / U*Vr~ l dxdv, 



(33) 



and demonstrated how this fact can be used in a per- 
turbation approach to slightly non-slender tori. Blaes 
(1985) used this theory to calculate growth-rate of the 
Papaloizou-Pringle instability and, more recently, the 
same approach has been employed in calculation of the 
eigenfrequency corrections to the epicyclic modes (Blaes 
et al, 2007). 

A small viscous damping of the oscillation modes can 
be treated in a similar way. In addition, in this paper we 
extend the perturbation approach of Blaes et al. (2007) 
to the case of a general angular momentum distribution. 
In appendix 1 we formulate first-order perturbation the- 
ory for the eigenvalue problem (29). Due to the viscosity 
term on the right-hand side of the equation (29), the eigen- 
frequency of the v-th. mode is slightly changed from the 
value er" that corresponds to the inviscid flow to 



= <r<°> 



r<l) 



(34) 



(29) The first-order correction is given by 



with 
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(35) 



with W„ ^ being the eigenfunction corresponding to the 
v-th mode of the inviscid torus. 

Since the oscillations of the inviscid tori are not 
damped, the growth-rate of the viscous torus is given by 

T^aflolm^ 11 ]- (36) 
The damping corresponds to the negative values of 7. 

5. Special cases 

5.1. The WKBJ limit 

In the limit of short wave-lengths of the perturbation, 
the eigenfunctions can be approximated using WKBJ 
ansatz 



W oc exp 



k T dx ~ 



kydy 



(37) 



where k x and k y are the horizontal and vertical compo- 
nents of the wavevector. Blaes et al. (2006) identified two 
distinct classes of modes, acoustic and inertial that obey 
the dispersion relations 



a 2 = Le and a 2 = R 2 '^ 



- e> ^ 

respectively. While the frequency of the acoustic modes 
grows with increasing k = (k 2 + k 2 ) 1 ^ 2 , that of the inertial 
modes is always smaller or at most equal R. 

After calculations briefly summarized in Appendix 2, 
we find that the approximate expression for the growth 
rates of the acoustic modes is 



IV 1 + 2)^°°' 



while that of the inertial modes is proportional to 
- 2 ^ attn 



7 , 
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(39) 



(40) 



(n+l)Xr 

where A* is a typical wavelength of the mode. Therefore, 
while the high-frequency acoustic modes are always 
damped by the viscosity, the low-frequency inertial modes 
become unstable when the a-viscosity is introduced. 

5.2. Keplerian angular momentum distribution 

In the limit of the Keplerian angular momentum distri- 
bution, the characteristic frequency of the inertial oscil- 
lations is equal to the radial epicyclic frequency, R = Cu r . 
The Lane-Embdcn function does not depend on x and is 
a function of y only, 

(41) 
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One set of the oscillation modes that depends linearly on x 
have eigenfunctions given by the Gegenbauer polynomials 
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Fig. 2. Damping rates of the lowest-order modes of constan- 
t-angular momentum tori. The solid and dashed lines corre- 
spond to the second-order and third order modes, respectively. 
The third-order modes are labeled according to their j and k 
numbers and x and y parity (see Blaes et al., 2006). 

where j > 1 denotes a number of vertical nodes (see Blaes 
et al., 2006). The corresponding eigenfrequencies are 



(43) 



The vertical epicyclic mode corresponds to j = 1. 
Applying our theory we find that the growth rate of the 
j-th mode is 

a\ (j + 2n)(j 
2, 



7j 



1 



1 



u> 2 aQr] 



(44) 



2(n+l) 

(see Appendix 3). Hence, all the modes arc damped in 
the Keplerian limit. In addition, for high-order modes 
with j ;g> 1 the damping rate agrees with that predicted 
by the WKBJ approximation. 

6. Results for the lowest-order modes 

Blaes et al. (2006) derived eigenfunctions and eigen- 
frequencies of the lowest-order modes of slender tori with 
arbitrary angular momentum distribution. The eigenfunc- 
tions of these modes are given by the polynomials of low- 
order in x and y. The most simple is the corotation mode 
corresponding to a = and a zero-order polynomial W = 1. 
It may become unstable when (3 > and in that case it cor- 
responds to the prinicipal mode of the Papaloizou-Pringle 
instability. The eigenfunctions linear in x and y character- 
ize two epicyclic modes - in slender tori they corresponds 
to a rigid motion in radial or vertical direction with asso- 
ciated epicyclic frequencies. Four second order modes (the 
X-mode, inertial mode, breathing mode and plus-mode) 
have eigenfunctions quadratic in x and y. Their velocity 
patterns have single node at the center of the torus. Blaes 
et al. (2006) show how to construct modes of arbitrary 
order by solving appropriate sets of linear equations. In 
this section, we calculate viscous growth-rates of the low- 
est order modes (up to the third order) using the equations 
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(35) and (36). For this purpose we use formulae for the 
eigenfunctions and eigenfrequencies derived by Blaes et al. 
(2006). 

6.1. The epicyclic modes 

The radial and vertical epicyclic modes correspond to 
the linear eigenfunctions, Wrad = x or W ver t — y. In 
the inviscid slender torus, they describe uniform displace- 
ments from the equilibrium positions with frequencies 
equal to epicyclic frequencies of a freely moving particle. 
Therefore, the shear of the corresponding velocity field 
vanishes and so may vanish the viscous damping. In fact, 
this is the case of the vertical epicyclic mode for which 
FWvert = and therefore 7 VC rt = 0. In the case of the 
radial epicyclic modes, however, we find that 



FW, 



1 



rad 



(uj 2 -R 2 )(4-k 2 )x 



and therefore the viscous growth-rate is 
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7rad 
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(A-R 2 )aCl . 



(45) 



(46) 



This is because during the radial oscillations, the viscous 
force does not vanish even though the shear of the per- 
turbed velocity field does. Because of variations in the 
kinematic viscosity coefficient rj (that arise due to a pres- 
sure variations), a nonzero azimuthal viscous force ap- 
pears, 



^visc K 



2/3r 2 



(4 — k 2 ) exp[— i(u)t — m<f>)] ■ 



(47) 



This force is in the anti-phase with the velocity perturba- 
tions 



■ exp[— \{ujt 



(48) 



2t(oj 2 - R 2 )f3r 

and therefore causes the damping of oscillations. Similar 
effect is absent in the case of the vertical oscillations be- 
cause the ^-component of the shear of the equilibrium 
velocity vanishes. 

Equation (46) shows that the growth-rate vanishes in 
the limit of the Keplerian angular momentum distribu- 
tion. This is because the amplitude of the pressure vari- 
ations decays when approaching that limit (note that 
Sp/p^ (1 — to 2 /R 2 )(Sv^ /^l) for the radial epicyclic mode). 
This is consistent with the behavior of p- modes in geomet- 
rically thin accretion disks. Although their amplitudes 
grow, the growth-rates are inversely proportional to the 
squared wavelengths of the modes (e.g. Kato, 1978; Kato 
et al., 1998). The p-modes with infinite wavelength (that 
corresponds to our radial epicyclic mode) are not affected 
by the viscosity. 

6.2. Constant-angular-momentum tori 

Figure 2 shows a radial dependence of damping rates 
(in units of afio) of the second-order (solid line) and third 
order modes (dashed lines). The third-order modes have 
been labeled according to the convention of Blaes et al. 
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Fig. 3. Damping rates of the plus mode and the breathing 
mode for different values of the bulk viscosity. The solid lines 
correspond to the case a = — 2/3 (£ = 0), the dashed lines are 
for a = —1/3,0 and 1/3 (£ = r)/3, 2r)/3 and r?, respectively). 

(2006). The coefficient a was set to a = —2/3, what corre- 
sponds to vanishing bulk viscosity coefficient, £ = 0. One 
may notice that the gravitational field enters into all for- 
mulae only through the expressions for the orbital and 
epicyclic frequencies fio, ui r and L) z . Similarly to other 
authors, we use the Schwarzschild expressions for these 
frequencies, despite the fact that the formulae have been 
derived in the framework of Newtonian physics. We be- 
lieve that this inconsistency will not affect the qualitative 
discussion of the paper significantly. The radial coordinate 
of the torus center is shown in units of the gravitational 
radii GM/c 2 . 

The analytic formulae for damping rates are rather 
complicated. The exception is the X-mode, whose 
growth/damping rate is 



2Qi 1 



7x 



■ aflo- 



(49) 



The damping rate does not depend on the bulk viscosity 
coefficient a because the X-mode oscillations are incom- 
pressible. On the other hand, both plus and breathing 
modes weakly depend on the value of a as it is shown in 
Figure 3. 

Figure 2 suggests that the the ratios j v /(aClo) tend to 
the constant for r ^> GM/c 2 . This behavior nicely illus- 
trate the scaling invariance of the Newtonian gravity. The 
Newtonian value of the damping rates of the plus and the 
breathing modes are 7+ = 2aQo and 7breath = (2 + a)af2o- 

6.3. Dependence on the angular-momentum profile 

The behavior of the damping rates of the four lowest 
order modes with changing angular momentum distribu- 
tion is shown in Figure 4. The top panel is devoted to 
the situation when torus is placed close to the central 
object at radius ro = lOGM/c 2 . The lower panel shows 
the same for the torus at large radii and corresponds 
to the Newtonian limit. The damping rates of all the 
modes except the breathing one are significantly reduced 
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Fig. 4. Damping rates of the lowest-order modes as func- 
tions of the angular momentum distribution. The top panel 
corresponds to the torus placed at ro = 10GM/c 2 , the bot- 
tom panel shows Newtonian limit r — > oo. In both situations 
we set a = — 2/3 what corresponds to vanishing bulk viscosity 
coefficient. 

for Kcplcrian angular-momentum distribution. Although 
our analysis predicts that the damping of all three modes 
vanishes completely, it is still possible that the actual 
damping rate is of a higher order in /?, say on the timescale 

of > (^any 1 . 

The growth-rate of the breathing mode in the Kcplcrian 
limit (ft —> 1) is 

7breath = ( X + ^) "^0- (50) 

Our results also suggest that the torus is unstable with 
respect to the inertial mode. In the relativistic case (top 
panel) this mode grows in time independently on k, while 
in the Newtonian limit this growth is suppressed for less 
steep angular momentum profiles. Figure 5 shows regions 
of stable and unstable configurations in the (r, ft)-plane. 

7. Discussion and conclusions 

In this paper we have studied viscosity effects in slen- 
der tori. Due to the angular- momentum transport, the 




20 40 60 80 100 



MGM/c 2 ] 

Fig. 5. Instability threshold for the inertial mode. 

torus becomes wider in the horizontal direction and thin- 
ner in the vertical as the angular momentum distribution 
approach the Keplerian one. We have shown that this 
process occurs on the timescale ~ a^ 1 ^ 1 and can be re- 
garded as a quasi-steady, i.e. as a series of steady inviscid 
tori configurations, because the poloidal velocity due to 
the viscosity is negligible. We have also found that the 
angular momentum transport is symmetric, the location 
of the center of the torus remains unchanged. This result 
is related to the reflection symmetry around the x = 
surface that appears in the slender torus limit. Including 
the higher-order terms in f3 that break this symmetry, we 
would likely find a shift of the torus center towards the 
central object. 

In the rest of the paper we dealt with the oscillations of 
the viscous tori. For simplicity, we neglected the secular 
evolution what allowed as to assume the harmonic time 
dependence of perturbations. We derived inhomogeneous 
master equation (29) with left-hand side giving the normal 
modes of the inviscid tori and right-hand side describing 
effects of viscosity on the modes. We formulated a gen- 
eral first-order perturbation theory for the case when the 
right-hand side is small and with the aid of this theory 
we obtained a general formula for the mode growth rates. 
Applying this formula in the short-wave limit, we have 
found that the viscosity acts to damp the acoustic modes. 
On the other-hand the viscous tori are unstable with re- 
spect to the inertial oscillations. The growth-rate of the 
instability increases with decreasing wavelength. We also 
calculated the damping rates of the acoustic oscillations 
in the Keplerian limit and showed that they agree with 
the results based on the WKB J theory in the limit of high 
nodal numbers. Finally we calculated the growth rates of 
the lowest order modes and explored their behavior with 
changing various parameters. We found that the damp- 
ing of three of four modes is significantly reduced as the 
angular momentum distribution approaches the Keplerian 
one. 
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Appendix 1. First-order perturbation theory for 
tori with arbitrary angular-momentum distribu- 
tion 

Here we deal with a perturbation theory for equation 
(29) with the assumption that a is small. The right-hand 
side is therefore treated as a perturbation. We start with 
the observation that both, the operators B and C are 
Hcrmitian with respect to the scalar product (33), i.e. 



(Al) 
(A2) 



The eigenfunctions that solve equation (29) with 

a = form a complete set and can be therefore used as 
a basis with respect to which solutions of the perturbed 
equation will be expanded. However, contrary to the 
first-order eigenvalue problems, the eigenfunctions Wi! 
are not orthogonal. Instead, they satisfy so the pseudo 
orthogonality relation 



(0) 



+ 



Wj?\BW^)=0 



(A3) 



for \i^v. 

The solutions of the perturbed problem are then ex- 
panded via a as 



aa 



(i) 



(A4) 
(A5) 



and substituted into equation (29). Comparing coefficient 
of the same powers of a we obtain in the first order 



(0)2 



.(0)2 



-M2 + -(0)2\ + £ 



w?, 



+2a^a^ [2a^ 2 +B] W^=iFW^. (A6) 

If we now perform scalar product with W„ and use the 
pseudo-orthogonality relation, we get 

2a^a^(w^\(B + 2a^)w^ 

i(w?\FWW), (A7) 

from which we imediatelly recover the relation (35). 

Appendix 2. The WKBJ approximation 

In the WKBJ limit, k — > oo, we find that 
(a 2 ~R 2 )f 2 k 2 



FW 
and 



4n(n4 
(2 + a)k 2 



I) 



(3 + a)S 2 



-2a 2 )W = 



-77- A 2 
2n 



2a z 



W. 



(A8) 



(A9) 



For the acoustic modes, the second term in the square 
bracket of equation (A8) vanishes because a 3> R and using 
the corresponding dispersion relation, we obtain 

-(W,W). (A10) 



W,FWj 



(2- 



1 



Similarly, we find that 



W,{B + 2a 2 )W) =a 2 (W,W) . 



(All) 

With aid of the equations (35) and (36), we therefore find 
that the growth/damping rate is 



it 



7 : 



1 



n+1 

In the case of the incrtial modes we find that 



W,FW) 



K 



-.2 



a 



4n(n + l)er 



and 



1 



W,{B + 2a 2 )W' / 
The growth rate of these modes is therefore 



--{W^fW). 



7 : 



(W,k 4 f 2 W) 



(A12) 



(A13) 



(A14) 



(A15) 



2(n+l)a 2 (W,k 2 fW) 

Both scalar products on the right-hand side are positive 
as their are integrals of the positive quantities. Therefore 
they cannot change the sign of 7 and we can conclude that 
the high-order inertial modes grow. Moreover, introduc- 
ing a typical wavelength of the mode A* , the upper and 
the lower scalar product can be estimated as ~ X~ A (W 1 W) 
and ~ A~ 2 (W, W) , respectively. Hence the growth rate can 
be estimated as 



No. 



1-4 
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7r(2n - l) 2 (2n + l) 2 ui 3 z T(j + 2n + 1) 



R 2 J (n+l)Aj' 



(A16) 



4"+ 1 (j + n- l/2)r 2 (n + 3/2)r(j - 1) ' 
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(A25) 



. j- o t-> • ^ j-u »i i i • . Then, substituting- into equations (A18) and (A20) and 
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Keplerian limit 

In the Keplerian limit, the Lanc-Embdcn function / is 
independent of x and the eigenfunctions Wj depend on 
x at most linearly. Applying the operator F given by 
relation (32) on W we obtain 

(a 2 -u> 2 r )(2 + a) d 2 



FW- = -■ 

3 Anin + l^jf^dy 2 



dy 2 



(A17) 



Performing twice the integration by parts in y, we find 
that the scalar product (Wj,FWj) is given by 



(Wj,FWj) 
where 



I 2 = 



d 2 W, 



dy 2 



(a 2 -^)(2 + q) j2 
An(n + l)a 



f n+L (y)dy, 



(A18) 



(A19) 



where we omitted the factor coming from the integration 
over x. Similarly, 

1 



W h {B + 2a*)W j ) = -—I 1 + (2a 2 - fi 2 )I° 



where 
7° = 
and 



\w 3 \ 2 p-\y)dy, 



f n (y)dy. 
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(A20) 
(A21) 

(A22) 



Using the substitution £ = uj z y and with the aid of stan- 
dard rules for derivatives of the Gegenbaucr polynomials 
and their normalization, the integrals 1°, I 1 and I 2 can 
be evaluated as 



1 

w 2 j 



1 - ,2 

-1/2 



C 



7iT(j + 2n-l) 



w z 4™-!(j + n+ l/2)r 2 (n - l/2)T(j + 1) : 
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then finally into the equations (35) and (36), we obtain 
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